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Abstract We explain how Regge theory and perturbative evolution may be made compatible at small 
X. The result not only gives striking support to the two-pomeron description of small-x behaviour, 
but gives a rather clean test of perturbative QCD itself. When x is very small, the proton's gluon 

CN . distribution function is significantly larger than is commonly believed. Perturbative evolution is invalid 

g; below Q2 ~ 5 GeV2. 
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1 Introduction 

There is a growing realisation 

[1][2][3] 

that the conventional approach to perturbative evolution breaks 
down at small x. Certainly, when, at a given fixed small x, the structure function varies rapidly with 
I Q^, it is not sufficient to expand the DGLAP splitting function in powers of and truncate the 
expansion after one or two terms. Rather, a resummation is needed and, while there are models of 
, how to do this, so far no reliable method is available. 

The Regge approach is apparently orthogonal to that of perturbative evolution. We have shownt^] 
that Regge theory gives a very good description of the data, not only for the proton structure function 
F2{x,Q'^), but also for its charm component F2{x,Q'^) and for J/tp photoproduction. The data 
' indicate rather clearly the need for a second pomeron, whose trajectory has intercept (1 + eo) ~ 1.4, in 
addition to the familiar soft pomeron with intercept (1 + ei) about 1.08. We call the second pomeron 
the hard pomeron. While at one time it was hoped that one might calculate its intercept from the 
BFKL equation, it now seems likely[5] that such a calculation is not within the scope of perturbative 
• QCD. As we explain in this paper, perturbative QCD merely governs how the magnitude of the hard 
pomeron's contribution to the structure function increases with Q^. 

X ■ Our belief ist^l W that the hard pomeron is already present in real-photon amplitudes and contributes 
. with the same fixed power {W'^Y" at all values of Q^. That is, at sufficiently large W 

i^2(a:,Q')~/o(g')x-'^° (1) 
for all values of and, again at sufficiently large W, 

a^P ~ 47r2aEM^o {W^Y" (2) 

with 

X, = iQ'r'-'^uQ') ^ (3) 

There is no theoretical reason why the behaviour should be a power rather than something more 
complicated, but we have foundW that a power fits the data very well. If Xq is to be finite, /o(Q^) 
must vanish as {Q'^Y'^"" at = 0. A very good fit to data, all the way from = to 5000 GeV^, 
is provided by the economical parametrisation 

/o(g2) = Xo (Q')^+^V(i + QVqIY^'^'" (4) 
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with Qo ~ 3 GeV. In this paper, we show that this form agrees remarkably well with what is obtained 
from DGLAP evolution, over a large range of . At present, it is not possible properly to apply 
perturbative evolution to the soft-pomeron component of F^ix^ Q^)- 



2 The DGLAP equation 

The singlet DGLAP equation 



^u(x, Q2) = I' dz P(z, a,(Q2)) u(|, 



where P is the splitting matrix, t = log(Q^/A^) and 



(5) 



Write the Mellin transforms ^ 

Jo 

P{N,a,{Q^))= [\zz''P{z,a,{Q^)) (7) 
Jo 

Then 

^u{N, Q2) = p(Ar, ^(^^ q2^ 

A power contribution (1) to F2{x,Q'^) corresponds to a pole 

in u(A^, Q^). More generally, consider a contribution x~'^°f{x, Q^). We assume that f{x, Q"^) vanishes 
at a; = 1 and is differentiable, for all Q^. Insert in the Mellin transform integral and integrate once by 
parts, to get 

-— ^ / (10) 

-fV — eo Jo 

So there is a pole at AT = eo with residue 

- / dx/,(ar,g2)=/(o,Q2) (11) 
Jo 



With 4 active quark flavours and a flavour-blind hard pomeron, fq{Q^) = ^fo{Q^)- We find[6] on 
taking the residue of the pole at AT = eo on each side of the Mellin transform (8) of the DGLAP 
equation, 

^f(Q2) = P(Ar = eo,a,(Q2))f(Q2) (12) 

We have previously!^] fitted accurate ZEUS and HI data in the range x < 0.001, 0.045 < < 35 
GeV^, together with data for a'^P. The result is 

eo = 0.437 (13) 
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Figure 1: ZEUS datal^] for the charm structure function F2{x,Q'^) and the charm photoproduction 
cross section, with the prediction from the fit to F2{x,Q'^) assuming that the hard-pomeron is 
flavour bhnd. 



and 



9.11 GeV^ 



Xq = 0.00146 



(14) 



The error on each of these quantities is large, but we have given their values to this accuracy because 
their errors are strongly correlated. 

If we include 4 flavours of quark and antiquark in the sum in (6), then at = 20 GeV^ the singlet 
quark distribution x^j^{qf + qf) ~ 0.095x~'^'* at sufficiently small x. We found also that the charmed- 
quark component F2 of F2 is apparently governed almost entirely by hard-pomeron exchange at small 
X, even at small values of Q'^, and that, within the experimental errors, its magnitude is consistent 
with the hard pomeron being flavour-blind. The hard-pomeron description of F2 is shown in figure 
1. According to perturbative QCD, the charmed quark originates from a gluon in the proton and the 
two distributions are proportional to each other to a good approximation over a wide range of x and 
Q^t^]. So this implies that the gluon distribution also is hard-pomeron dominated. The conventional 
approach to QCD evolution is correct if x is not too small, because it does not then probe the splitting 
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matrix at z = 0, which is where its expansion in powers of is illegal. We assume that at x = 0.01 
and = 20 GeV^ the value of g{x,Q^) extracted by the HERA experiments is reasonably close to 
the correct value. At Q'^=20 GeV^ and x = 0.01, an NLO fitl^O] to the combined ZEUS and HI data 
gives xg{x,Q'^) = 5.7 ± 0.7. Other authorsl-'^-'^]!-'^^^!^] gj^^j much the same value. This is 8 ± 1 times 
hard-pomeron component of the singlet quark distribution. 

An unresummed perturbation expansion of the splitting matrix P{N,as) is not validt^l for small 
values of N because it introduces a spurious singularity at = 0. According to (12) and (13), we 
need P(A, a^) at a value of N far from 0, and so it is reasonable to hope that resummation is not 
needed. The numerical values of the elements of the matrix P(N,as) in one and two-loop order are 
plotted in figure 4.9 of the book by Ellis, Stirling and Webbert'^] for the value = 27r/30. Prom these 
we may evaluate the running-coupling splitting matrix P[N,as{Q^)) in one and two- loop order*. We 
use the one and two-loop forms for the running coupling!^] 

where 

b=''-'^f 6' = ^^^^^ (156) 
127r 27r(33-2n/) ^ ' 

In each case, we choose A such that as{M^) = 0.116. This gives 

A^*^ = 140 MeV A^^*^ = 400 MeV (16) 



We use four flavours throughout as, at the energies we are considering, the charm contribution is 
active and the beauty contribution is so small that its omission has a negligible effect. 

We integrate the differential equation (12). The result for the singlet quark distribution is shown 
in figure 2a, where the solid curve is the result of the two-loop-order perturbative QCD evolution 
according to (12), and the broken curve is the fit to the data given by (4). We have taken the ratio of 
the gluon distribution to the hard-pomeron component of the singlet quark distribution to be 8.0 at 
= 20 GeV^. Note that the gluon/quark ratio is a parameter which in principle we could change. 
However it turns out that the value of 8 we have obtained from the NLO fit to the combined ZEUS 
and HI data works well. Figure 2b shows how the gluon distribution 

xg{x,Q^) = fg{Q^)x-'° (17) 



evolves. Figure 3 shows that there is very little difference between one-loop-order and two-loop-order 
evolution, except at small Q^. This is because we have used the splitting function for a value of N 
safely away from N = and it encourages the hope that resummation, if we knew how to perform it, 
would make little difference to these results. 



3 Discussion 

We have a number of comments on these results: 

1 It is evident from figure 2a that two-loop perturbative QCD describes the evolution of the strength 
of the hard-pomeron contribution to F2{x,Q^) extremely well for greater than about 5 GeV^, up 
to values of Q"^ beyond where data exist. The perturbative QCD differential equation is a large-Q^ 



* It is not made clear in the book that the qg plot includes the necessary factor 2n / . 
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Figure 2: (a) NLO evolution with A = 400 MeV of tlic hard-pomcron coefficient /o(Q^) (solid curve 
and the fit (4) (broken curve); (b) evolution of the gluon structure function xg{x,Q^) at x = 10~^. 
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Figure 3: NLO evolution with A = 400 MeV (solid lines) and LO evolution (broken lines) with A = 140 
McV (a) of the hard-pomeron coefficient fo{Q^) and (b) of the evolution of the gluon structure function 
at X = 10-5. 



equation and is not valid when is small; it even makes the hard-pomeron contribution to F2 become 
negative below 2.5 GeV^. 

2 Our fit to the data, from which we extracted the hard-pomeron contribution, used data from Q"^ = 
to 35 GeV^. Over most of this range, the hard-pomeron contribution is only a small part of F2{x, Q^), 
though its relative magnitude has increased by = 35 GeV^: see figure 4. So, although it is 
well-establishedl-'^-'^'-'^^] that the complete F2{x,Q'^) obeys perturbative QCD evolution beyond = 5 
GeV^ for the range of x where data exist, it is not at all trivial that the hard-pomeron part of it does. 
This successful link with perturbative QCD provides rather striking verification of the hard-pomeron 
concept and of our having extracted it correctly from the data. 

3 It also provides a clean test of perturbative QCD itself, one that genuinely tests evolution rather 
than being a global 

fit[ll'12]. 

It depends on just one parameter, the ratio of the gluon distribution to 
the hard-pomeron part of the singlet quark distribution at some value of Q^, which we chose to be 
20 GeV^. The curve in figure 2a uses the value 8.0 for this ratio but, as we have said, analysis of the 
HERA data gives an error ±1 on this number. Changing it within this range has a fairly small effect 
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Figure 4: The thick curves are fits to F2{x,Q'^) at = 5 and 35 GeV^, with the HI data used at 
those values of Q^. The thin curves are the hard-pomeron contributions. 



at large Q^, about 10% at 5000 GcV^, but causes quite a large change in the details of how it breaks 
away from the phenomenological curve when we evolve back to small Q^. Similar remarks apply to 
changing A within its allowed range of values. 

4 Our phenomenological fit behaves at large as the power , while the solution to the perturbative 
QCD differential equation (12) rather behaves asymptotically as a power of logQ^, approximately 
(logQ^)^; nevertheless the two agree well over a very large range of Q^. According to figure 2a, they 
begin to diverge from each other only at values of beyond the data, so in our fit we could instead 
use the perturbative QCD form. However, it is very much more complicated to parametrise and, more 
importantly, using the perturbative QCD form at small is incorrect. Using different forms for 
different ranges would lead to matching problems and it would be almost impossible to achieve 
the required analyticity in at the join. 

5 The conventional approach to evolution expands the splitting matrix P{N,as) in powers of as- 
Successive terms in the expansion are increasingly singular aX N = 0. This is a signal that the 
expansion is illegalt^J for small values of A'', since the complete splitting matrix is regular at A^ = 0. 
Because an unresummed expansion that needs the splitting matrix at small A^ makes the splitting 
function larger than it really is, a gluon distribution of a given magnitude apparently gives stronger 
evolution than it really should. That is, the conventional approach will tend to under-estimate the 
magnitude of xg{x,Q'^) in certain regions of {x,Q'^) space. This is verified by our results for the 
evolution of xg{x,Q^): figure 5 shows the proton's gluon structure function at two values of Q^, 
according to our calculations, which do not use the splitting matrix at small A", and compares it with 
what is extracted from the data by conventional means. 

6 We cannot use a similar approach to the evolution of the contribution from the soft pomeron, 
because this needs P{N,as) at N ^ 0.08, dangerously close to A^ = 0. Handling this will need 
resummationt"'^^] ["^1 1"'] , so as to tame the singularity at A^ = 0, but we do not yet have a reliable way to 
perform this resummation. As we have explained, experiment finds that the charm-quark distribution 
does not contain a soft-pomeron term; we are assuming that our inability to handle the soft pomeron 
is therefore not relevant for the gluon distribution. 

7 The fact that we can fit the data so well with the only singularities of u{N, Q^) in A" > identified 

with the hard and soft pomerons, supports our behef[6] that P ( A^, as ) does not have any singularities 
in this region. This again contrasts with the conventional attitude: we maintain that the relevant 
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Figure 5: Gluon structure function xg{x,Q'^) at (a) 



(b) 

20 and (b) 200 GeV^. 



0.001 



In each case the 



thick hue is our evolved distribution. In (a) the thin hnes are the hmits extracted by conventional 
NLO analysis of HERA datal^^^l In (b) the middle line ist^^.lS] CTEQ5M and the lower line ist^^'l^l 
MRST20011. 

singularities of u(A^, Q'^) in the complex A'^-plane are not generated by the perturbative evolution, 
which merely governs how the strength of their contribution increases with Q^. They are already 
present at = and their position does not vary with Q^. We have explained before!^] that there 
are strong theoretical reasons to believe this: one cannot simply start perturbative QCD evolution at 
some value of and ignore how this joins on to what is found for smaller values. 

8 We stress again how very few parameters we need to fit the data for F2{x,Q^) for x < 0.001. The 
hard-pomeron contribution is parametrised in terms of the three parameters Xq, cq and Qq introduced 
in (4), and for the soft pomeron we have similarly^ Xi,ei and Q^: 



(18) 



Of these, ei 0.08 is determined from soft hadronic interactions, and Xi from the fit to a'^'^. The 
latter requires t^l also a contribution from /2 and 02 exchange, but this does not contribute very much 
to F2{x, Q'^) for X < 0.001. 

9 We have remarked that, although the data for F2{x,Q'^) are now highly accurate, the value of eo 
is still quite poorly determined. It depends!^] on what is assumed for the large-Q^ behaviour of the 
soft-pomeron coefficient function fi{Q'^). In this paper we have used eo = 0.437, which corresponds 
to asuming that /i(Q^) has the form (18) and so goes to a constant at large Q^. But one may obtain 
an equally good fit to the data for F2{x,Q'^) by assuming that /i(Q^) vanishes like 1/Q at large 
Q^. Figure 6a shows the equivalent of figure 2a for this case, with eo = 0.394, and figure 6b is the 
equivalent of figure 5a, from which it is seen that the gluon structure function needs to be slightly 
larger at = 20 GcV^ than in the previous case. By = 2000 GeV^ it is quite a lot larger: see 
figure 7. The lower curve in this figure is the CTEQ4M prediction, based on conventional unresummed 
NLO evolutionl2. 

10 The LO or NLO perturbative QCD evolution makes the coefficient fo{Q^) of the leading power of 
1/x increase indefinitely with increasing Q^. In order not to conflict with the momentum sum rule, 

at any fixed x nonleading powers of x must become progressively more important, so that the largest 
value of x for which the leading power alone gives a good approximation to F2{x,Q^) decreases as 
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Figure 6: Equivalents of (a) figure 2a and (b) figure 5a for the case eo = 0.394 
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r igure 7: Gluon distribution at = 2OOO GeV^. The upper line is the perturbative QCD-evolved 
distribution for eo = 0.394, the middle line for eo = 0.437, and the lower line is the MRST20011 
two-loop predictionl^^'-'^'^J 

increases^. This is why, while our plots of the gluon distribution at = 20 GeV^ extend up to 
X = 0.01, for 200 GeV^ we stop at x=0.001, and for 2000 GeV^ at 0.0001. We guess that these values 
for the limits of the validity of the single-term approximation are safe. 

11 To summarise, what we have achieved is a genuine evolution, in contrast to global 

fits[ll'12] 

which use data at large to help constrain the fit already at small Q^. Wc used data up to 
= 35 GeV^ to determine the hard-pomeron component of F2(x, Q^); in this fit we did not assume 
perturbative QCD evolution, but just made a numerical fit to the dataW. Now we find that, for the 
hard-pomeron contribution, we have almost perfect agreement with perturbative QCD evolution up 
to = 5000 GeV^. This is a striking verification that the hard-pomeron concept is correct, as well 
as being a success for perturbative QCD itself. It leads us to conclude that, at very small x, the gluon 
structure function is somewhat larger than has until now been believed. 
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